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Generic quantum systems -as much as their classical counterparts- pass arbitrarily close to their initial state 
after sufficiently long time. Here we provide an essentially exact computation of such recurrence times for 
generic non-integrable quantum models. The result is a universal function which depends on just two param¬ 
eters, an energy scale and the effective dimension of the system. As a by-product we prove that the density of 
orthogonalization times is zero if at least nine levels are populated and connections with the quantum speed limit 
are discussed. We also extend our results to integrable, quasi-free fermions. For generic systems the recurrence 
time is generally doubly exponential in the system volume whereas for the integrable case the dependence is 
only exponential. The recurrence time can be decreased by several orders of magnitude by performing a small 
quench close to a quantum critical point. This setup may lead to the experimental observation of such /a5t 
recurrences. 


I. INTRODUCTION 


The idea that the world will repeat itself after some time 
is an ancient one and appears in many philosophies and cul¬ 
tures. Babylonians called the Great Year the time needed for 
the planets to return to their initial conditions iniEi El. How¬ 
ever it was only in the end of the XIX century that Henri 
Poincare rigorously proved the existence of a recurrence time 
for a certain description of reality a. The beauty and the 
power of Poincare’s result lie in its simplicity and in its wide 
range of applicability. The only hypotheses of Poincare re¬ 
currence theorem are that 1) the (phase) space available to the 
dynamics is bounded in volume; and 2) the dynamical flow 
p* preserves the volumes 0. Under such hypotheses the sys¬ 
tem returns arbitrarily close to its initial conditions after suf¬ 
ficiently long time. Hypothesis 2) is satisfied if we believe a 
classical, Hamiltonian, description of the world thanks to Li- 
ouville’s theorem. Hypothesis 1) holds true if, for instance, 
the motion is restricted to a bounded region of space because 
of conservation of energy 0. Poincare recurrence theorem 
has philosophical and counterintuitive implications. For ex¬ 
ample -if we believe a conservative, classical dynamics- the 
theorem predicts that, if we let two different gases mix by re¬ 
moving a barrier, there will exist a time where we will see 
the two gases separate again. Indeed Zermelo used Poincare 
result’s to argue that Boltzmann’s formula for the entropy 
should decrease after some sufficiently long time, thus vio¬ 
lating the claim that entropy always increases Q. The reason 
why we never see the two gases nicely separate back again is 
due to the astronomically large times needed to observe such 
recurrences. Poincare’s argument does not give any indica¬ 
tion on how to estimate such recurrence times. In order to do 
that one must generally have stronger assumptions. The very 
first estimate of recurrence time has been given by Boltzmann 
himself in reply to Zermelo’s criticism 0. Boltzmann esti¬ 
mated that the time needed for a cm^ of gas to go back to its 
initial state, has of the order of many trillion of digits. 

In more recent years the estimation of such Poincare recur¬ 
rences became an active research topic belonging to the field 


of dynamical systems 0. The first estimation of return times 
appeared in lITOl for a system of harmonic oscillators. More 
precisely in Ref. ifTOll the average recurrence has been com¬ 
puted. It is defined by Tp = lim,r->oo tr/r, where tr is the 
location of the r-th recurrence. As noted by Kac in HU, the 
result of Ref. HO) can be obtained almost immediately thanks 
to the following formula of Smoluchowski, valid for a discrete 
evolution, 


H{S) - n{A) 

^ A)' 


( 1 ) 


The above formula gives the average time taken to return in set 
A given that we started from a point of the set. In Eq. ([^ p{S) 
is the (finite) invariant measure of the total available phase 
space and g'^ is the mapping advancing one unit of time t (see 
e.g. US). The result for continuous evolution is obtained tak¬ 
ing the limit r —> 0. Eq. Q can readily be applied to a clas¬ 
sical integrable system. The assumption of finiteness of the 
available phase space (more precisely the manifold is assumed 
to be compact and simply connected), implies that the motion 
takes place on a dimensional torus = {(^i,..., (j)N) 
mod 27r}, where N is the number of degrees of freedom ifTSl . 
The map is given by g^tp = {pi -F wif,..., (j)iq -\- 
If the frequencies {wi} are rationally independent (RI, i.e. lin¬ 
early independent over the field of the rationals), time aver¬ 
ages are equivalent to phase space averages on the torus and 
Eq. Q applies. In this case, if A is taken to be the angular in¬ 
terval {Api ,..., Ap^), Eq. 0 gives, after taking the r —>■ 0 
limit, 


Tr = 


(nf., Sk - 0 

E N ujj 

j^l A^j 


( 2 ) 


Eor example, taking for simplicity Apj = ewe obtain Tr ~ 

We know however reality is ultimately quantum. Does a 
recurrence theorem applies also to the quantum world? That 
this is the case was shown in Ref. Da but no actual estimate 
of recurrence time was given. Some approximate estimates 
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for some specific initial states where given in lITSlfThl . In this 
paper we solve this problem in an essentially exact way. 

Let’s assume then, from here on, that the world evolves ac¬ 
cording to Schrodinger equation and that, for simplicity, the 
initial state is pure. For the time being we assume that the 
system’s Hamiltonian H has purely discrete spectrum. This 
is a crucial assumption on which we will come back (and 
weaken it) later. Then H = ErJln with orthogonal pro¬ 
jectors n„ and let !■!/)) = |'!/’(0)) the state of the system at 
t = 0. We define |n) = n„|'!/>)/ ||n„|'(/i)||. The state at time 
t is then {-ipit)) = = J2n l|nn|^/’)|| Ull. 

As time goes by {ipit)) spans a d dimensional torus (d, pos¬ 
sibly infinite, is the number of non-zero ||n„|'(/i)||) with radii 
11 n„ I'!/))! I and angular frequencies En- If the energies En are 
rationally independent (which is natural to assume) the torus 
is filled uniformly as time goes by. It is tempting then to try 
to recycle the classical result and use Eq. Q to estimate the 
recurrence time in the quantum case. The problem with that 
is that the distance of points on the torus does not imply any 
physical distance between states. For example states \n) with 
large amplitudes ||n„|'!/>)|| should count more than states with 
small amplitudes. A natural distance between {tpit)) and 
instead, is encoded in the fidelity 

Eit) = Mm)f- (3) 

The quantity in Eq. Q is also known under different names 
such as survival probabilities or Loschmidt echo, and it has 
been studied in various contexts such as quantum chaos, the 
theory of fermi edge singularities, equilibration and decoher¬ 
ence EUSl. The statistical properties of E{t) when t is 
seen as a uniform random variable on the infinite half-line 
have been investigated in a series of works lim 12311271. If 
the spectrum of H is discrete, E{t) is an almost periodic 
function and E{t) evolves indefinitely without ever con¬ 
verging to a limit as t —>■ oo. Conversely if the spectrum is 
purely continuous (or if only this part of the spectrum is popu¬ 
lated) limt_>oo E{t) =0 essentially because of the Riemann- 
Lebesgue lemma. In this situation the system does not return 
close to its initial state and actually becomes orthogonal to it 
for long times Il29l . Hence only if there is a discrete spec¬ 
trum, and if the corresponding probabilities = ||n„|'(/))|l^ 
are non-zero, the system has a chance to return close to its ini¬ 
tial state. We henceforth set ourself in this situation. In this 
case, typically, E{t) start from its maximum one at f = 0, 
drops to an average value in a time that may well be called 
equilibration time f23[ l30l . and starts oscillating around this 
value (see Fig.[T]for a typical plot). We are interested in esti¬ 
mating the time E{t) will go back to a value u G [0,1] possi¬ 
bly (but not necessarily) close to one. Let us proceed as fol¬ 
lows. We call tn the n-th solution of the equation E{t) = u 
on the half-line f > 0 starting from the left. We define the 
(average) recurrence time as Tfj(it) = lim„_>.oo This 

number clearly does not need to coincide with, for instance 
the first recurrence, but it’s essentially the only quantity that 
makes physical sense in that -as we will show- depends only 
on a finite number of physically relevant constants. Estimat- 
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Figure 1. Left panel: typical behavior of the fidelity for a many body 
system. The system is given by Hamiltonian Eq. {13} with K — 0 
with L = 50. The horizontal line is the location of the average 
The peaks observed at times approximately integer multiple of 
the system size are revivals discussed in im. Right panel: univer¬ 
sal function U{x) = (y7r/2)e“’/vT;. U{x) has essentially three 
regimes: for a; ^ l,h({x) ~ 1/^, for x = 0{l),h({x) ~ 0(1) 
whereas for x 3> l,U{x) ~ (dashed line). 


ing the exact recurrences f„, besides being mathematically al¬ 
most intractable, is physically meaningless. In fact the exact 
recurrences depend crucially on the energy levels En which 
must be known with infinite accuracy EiHia. Since the lat¬ 
ter seems a rather unphysical presumption we focus on the 
computation of Tn(u). 

Our strategy will be the following. First we estimate the 
number of zeros Nt(u) of the equation E{t) = it in a large 
segment [0,r] with T ^ 1. We will show that, for large 
T, Nt{u) ~ T D{u) with a finite density D(u). Then the 
recurrence time is given by Tfj(it) = limT^ooT/Nriu) = 
1 / 0 ( 11 ). 

We now show how to compute D{u). Similar techniques 
have been pioneered by Kac and Rice to study the average 
number of roots of a random polynomial or time signal ll^ 

Ea. 


II. NON-INTEGRABLE CASE 

We assume here the the Hamiltonian H has pure discrete 
spectrum and will consider generalizations later. This condi¬ 
tion plays an analogous role as the requirement of a compact 
phase space for the classical case. The Hamiltonian has then 
the following spectral decomposition H = with 

eigen-projectors n„. 

The number of zeros of the function F{t) — u in the interval 
[0, T], is given by the following formula 

Nt{u) = r 6{E{t) - u) \E'{t)\ dt, (4) 

Jo 

as can be checked by a change of variable (and with 
= dT{t)/dt). For large T we can use Nt(u) = 
TD{u) -\- 0{T^) where the density is given by D{u) = 
limT->.oo T~^ /o^ ~ u) \E'{t) \ dt provided D{u) is fi¬ 

nite. The density can then be computed with 

D{u)= I df\r\P^,^,{u,f) (5) 
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where /') := - f)S{J^'{t) - /') is the joint 

probability distribution of T and and overline indicates in¬ 
finite time average, i.e. a(t) := limT->.oo T~^ a(t)dt. For 
the next step we assume that the energies En are rationally 
independent, i.e. linearly independent over the field of ratio- 
nals. This is a very reasonable assumptions for generic non- 
integrable system. Thanks to rational independence we can 
express time averages as phase space averages over a large 
multi-dimensional torus US). The steps are detailed in the 
Appendix The main ingredient is given by the fact a cer¬ 
tain distribution related to Pjr jr,, becomes Gaussian in the 
limit of large dimensionality. 

The final result for the recurrence time is surprisingly sim¬ 
ple 




( 6 ) 


where E = ’YhjP] is the variance of the en¬ 

ergy {En} with respect to the ensemble Vn = Pn/'l2nPn’ 
i.e. The factor 

1//S.E is essentially needed to reproduce the correct dimen¬ 
sionality whereas the main dependence is through a universal 
function of x = u/E, U{x) = (y^/2)e®/-yx (see Fig. 
right panel). Since, forgetting the size dependence, Ai? ^ J 
where J is some energy scale of E[, a good approximation of 

Eq. ( 61 is given by Tr ~ J~^^ E . Several consider¬ 
ations can be drawn on the basis of Eq. 

1 ) Eor macroscopic systems the recurrence time is, in gen¬ 

eral, doubly exponential in the system’s volume. Indeed, in the 
typical situation, the average fidelity is exponentially small in 
the system volume, E ^ with a positive constant (see 

e.g. ESjl l |l36l. Sufficient conditions are i) the Hamiltonian 
is local and extensive, meaning El can be written as El = 
Yx with h{x) locally supported operators; together 

with ii) the initial state is exponentially clustering, i.e. the state 
(•) satisfies {h{x)h{y)) — {h{x)){h{y)) = O 

Hence, under these very general conditions, one has readily 
Tr{u) ~ (forrt > E). 

2 ) It is tempting to relate the recurrence time to the con¬ 

cept of entropy. With the aid of the equilibrium state p := 
\tjj(t)){tjj(t)\ one can write El 2 {p) = — InJ^, the 2-Renyi 
entropy of p. Recall the entropy inequality El 2 {p) < S{p) 
where S is the Von Neumann entropy of p, S{p) — — trplnp. 
The inequality becomes tight in the limit of totally mixed 
p, i.e. when E is small. Since, as we have shown, typi¬ 
cally E ^ e~°‘^ a possible proxy to the recurrence time (for 
M ^ is given by (omitting pre-factors) Tr ^ , a 

formula sometimes found in the literature (see e.g.. Ref. llJTll '). 

3) Equation (j^ and its inverse D{u) = \/Tr{u) should be 
handled with care at the extrema m = 0 or m = 1. Eor example 
Eq.^ predicts a finite density of exact recurrences (u = 1) 
but a zero density of orthogonalization times (i.e. times for 
which E{t) = 0). This is quite surprising given the fact that 
typically E{t) spends much of the time close to its average 


E which in turn is very close to zero. One can ask wether 
these results are correct or are a consequence of the Gaussian 
approximation used to derive (|^. The tiny density of exact 
recurrences (for u = 1) predicted by Eq. (|^, is indeed a con¬ 
sequence of the Gaussian approximation. In fact, for ratio¬ 
nally independent energies the only solution of E{t) = 1 is 
at f = 0 (see Appendix [C| for a proof). Instead, using this ap¬ 
proach without the Gaussian approximation we have been able 
to prove that indeed the density of orthogonalization times is 
zero (see Appendix [^, provided at least nine energy levels 
are populated and the energies are RI. 

4 ) Eor u close to 1 the roots of E{t) = u are very close 
to a maximum of E{t), hence they must come in close pairs. 
This means that essentially, for it « 1 a more precise estimate 
of the recurrence time is given by 2Tr{u). More importantly, 
on the basis of numerical simulations, we observed that, for u 
close to one, the solutions of E{t) = u come in very narrow 
spikes. In other words, even if a time t exists for which the 
system returns close to its initial state, at a later time t + At, 
with At very small, the system is likely to differ strongly form 
the initial state, and be nearly orthogonal to it if E is small. 

5) Eq. has been derived assuming a finite number d of 
non-zero pj. Nothing really dangerous happens taking the 
limit d —t oo, provided the spectrum stays discrete and nu¬ 
merable and AE and E stay finite and non-zero (E is finite 
because E < (V'lV') = !)■ In f^ct also in this case E{t) is an 
almost periodic function. In a more general situation H can 
also have a continuous spectrum. The characteristic function 
x{t) = ('01'0(f)) can then be written as 

X{i) + [ g{e)e~^^^de, (7) 

n •' 

where the positive function g{e) (this is the continuous part 
of the projected density of states, PDOS, {4’\5{H — e)|0)) 
has support on the continuous spectrum. If g{e)de is ab¬ 
solutely continuous, the Riemann-Lebesgue lemma guaran¬ 
tees that limt_>oo / 5 (e)e“*‘'^de = 0 . Hence, in this case, 
the second term of Eq. (|7]i contributes at most a finite num¬ 
ber of zeroes. A similar argument with small modification 
works for the fidelity |x(f)|^- Hence we deduce that Eq. 
works also in this more general setting with the caveat that 
the pnS refer only to the discrete spectrum and one may not 
have YnPn — 1- Note that Eq. ^ should be utilized with 
u G (0,M) with M = YnPn ^ I- Moreover the Gaussian 
approximation employed may become worse if the missing 
term / g{e)de becomes a sensible fraction of 1. 

The above results inspire some connections for the debate 
around the quantum speed limit if we consider Eq. (|^ for 
M ^ 1 ll38l - l4TI . In some situations (for example to esti¬ 
mate the computational power of the ultimate laptop or the 
universe one is interested in counting the number of 

total distinct (i.e., orthogonal) states that a “system can pass 
through in a given period of time” ll^ . Calling Tqsl the 
minimum time it takes for a system to go to an orthogonal 
state, it is known that Tqsl A (^/ 2 )/(( 77 ) — Eq). Ac¬ 
cordingly, in Ref. ll^ . it was estimated that the maximum 
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number of states, -/Vmax, a system can pass through in a time 
interval T, is given by -/Vmax = (2/7r)((_ff) — Eq)T. How¬ 
ever, in this situation where the initial state is specified, there 
is no need to compute A^max^ but rather Nt{u) is the num¬ 
ber we are after. Now, if T is sufficiently large, the number 
of times the state returns at fidelity u, is precisely given by 
Nt{u) ~ T/Tr = T/S.ElU{ulF). Some considerations are 
in order, i) Exactly at u = 0, as discussed in 3), one has 
actually Nt{u) = o{T) (i.e. lim 7 ’_>,oo N{u)/T = 0) for R1 
energies when at least nine levels are populated; ii) The es¬ 
timate A^max of Ref. 1^ qualitatively agrees with our N{u) 
in the region where exponential and square root terms are in¬ 
effective, i.e. for u/F ^ 1. In fact in this region we have 
Nt{u) ~ TAE. Note that u/F ^ 1 is likely the region 
where one would want to use these results, since, as we have 
seen, typically F is very small. However the relevant energy 
scale is not {{H) — Eq), (and not even the standard devia¬ 
tion y/ (iT^) — (Tf)^), but rather the standard deviation AE 
computed with the “squared” distribution ■ 


III. INTEGRABLE CASE 


The result Eq. ^ is valid assuming rational independence 
of the many-body spectrum. This condition is massively vi¬ 
olated for quasi-free integrable models whose Hamiltonian is 
quadratic in creation and annihilation Eermi/Bose operator. It 
is an interesting question per se to compare typical timescale 
of integrable systems with non-integrable ones 1441 . We con¬ 
sider hence a quasi-free system of fermions where the Hamil¬ 
tonian El is bilinear in Ci, (cj Eermi annihilation operators) 
and the initial state is Gaussian (i.e. satisfies Wick’s theorem). 
Eor simplicity we consider a one-dimensional geometry and 
assume that the fidelity can be written in the form 

F{t) =Y\_[^ - akSin'^(tek/2)] , (8) 

k 


sian and using Eq. Q one obtains 

Tr = exp [(Inu - In J')V(2 ct|)] , 
az' 

with coefficients given by 

In ^ 

k 

k 

4 ' = 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


and 'Zk, are smooth, bounded, function of ak,ek- 

Some remarks are in order. 

1) The sum over k in Eqns. (10 1 -( 12 1 , runs over an extensive 
number of terms. Accordingly (since the summand functions 
are bounded) we conclude that both In F and are extensive 
quantity. This implies that the recurrence time for quasi-free 
fermions is only exponentially large in the system volume. 

2) The behavior at the border m —0,1 must be handled 
with particular care. Similarly as in the non-integrable case, 
Eq. (|^ predicts a finite density of exact recurrences times and 
a zero density of orthogonalization times. The small, finite 
density of exact recurrences predicted by Eq. (j^ at m = 1, is a 
consequence of the Gaussian approximation. In fact a similar 
argument as for the non-integrable case shows that, if all 
are positive, F{t) = 1 only at f = 0. 

3) The exponential vs. the double exponential dependence 
of the recurrence time with the system’s volume may suggest 
that Tr may serve as a detector of integrability. However, 
any (no matter how small) non-integrable perturbation of an 
integrable model will result (with probability one) in a ratio¬ 
nally independent many-body spectrum and consequently a 
Tr given by Eq. (§1. As a consequence Tr will be discon¬ 
tinuous at the integrable point. A similar behavior has been 
observed for the temporal variances 


where {e^} is the one-particle spectrum of H, k is a quasi¬ 
momentum label, i.e. k = Tr{2n + 1)/L, withn = 0,1,... ,L, 
L is the number of sites and ak S [0,1]. This is the case for 
example for quenches of the one dimensional XY model, but 
can be valid more generally provided the Hamiltonian and the 
state are translationally invariant. In particular generalization 
of Eq. ([^ to Z?-dimensions is straightforward promoting k to 
a D-dimensional vector in the Brillouin zone. 

This is one of the rare cases where the calculation is harder 
for integrable systems because the need of passing to the one 
particle space introduces additional complications. Eor quasi- 
free systems it is known that the logarithm of the fidelity (in¬ 
stead of the fidelity itself) becomes Gaussian distributed is. 
It is then natural to consider the variable Z{t) = log F{t). 
We will then estimate the number of zeroes of the equation 
logD(f) = In It. Proceeding as previously we are now led 
to consider the joint probability distribution Pz^z'{z,z') — 
S{Z(t) — z)S{Z'{t) — z'). Eor large L, Pz,z' becomes Gaus¬ 


IV. EAST RECURRENCES 

The double exponential growth of the recurrence time with 
the system’s volume, poses serious questions on the possibil¬ 
ity of observing recurrence times in practical situations. On 
the other hand in some cases it may be possible to prepare 
the initial state such that p„ = 0 but for a few terms. Rabi 
oscillations may be seen as a limiting case of this situation. 
Another possibility is provided by a small quench experiment. 
In a quench experiment the system is initialized in the ground 
state of the system’s Hamiltonian with certain parameters Ai. 
The parameters are then suddenly changed and the system 
is evolved with parameters A 2 = Ai -f (5A. Whereas a gen¬ 
eral perturbation results in roughly as many excitations as the 
Hilbert’s space dimension, for sufficiently small quench am¬ 
plitude (5A, the number of excited quasi-particles is propor¬ 
tional to the system’s volume. A sufficient condition for small 
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Figure 2. Recurrence time in small quench experiments near crit¬ 
icality. The system is initialized in the ground state of Hamilto¬ 
nian CD with parameter ki , hi and then let evolve with the same 
Hamiltonian with parameters K 2 , /i 2 . Left panel: non-integrable 
case with ki = K 2 = 0.4. Other parameters are L — 12, , h 2 = 
hi -F 5h, 5h — 0.03, u = 0.98. The recurrence time is given 
by Eq. The vertical line signals the estimated critical point 
he = 0.218. Right panel: integrable case with ki = K 2 = 0. Other 
parameters are L = 12, ,h 2 = hi+ Sh, Sh — 0.03, u = 0.98. 
The recurrence time is given by Eq. 12} The vertical line signals the 
critical point he = 1. 

quench is given by |i5A| ^ where L is the linear sys¬ 

tem’s size and the D spatial dimensionality. As a consequence 
the system evolves effectively in a small Hilbert space whose 
dimension is roughly given by 1/J^ ll46l and the recurrence 
time can be diminished by several order of magnitudes. In 
fact one has ~ 1 — (c positive constant 0(1)) and 

hence Tfj ~ . However a further reduction is possible. 

If the system is close to a critical point, the effective dimen¬ 
sion 1 / is brought down to 2 — 3 and Tr becomes roughly 
,/“^0(l). A sufficient condition for small quench becomes in 
this case |(5A| <C where ly is the critical exponent of the 

correlation length ^ ED, whereas the quasi-critical regime is 
defined by by ^(Ai) L. In this regime the recurrence time 
is roughly size independent and may well become observable. 

To illustrate these effects we show numerical results for 
the one-dimensional Ising model in transverse field with addi¬ 
tional next nearest neighbor interaction. The, so called, TAM 
Hamiltonian, is given by 

L 

H = -Y^ + hal) , (13) 

i=l 

with periodic boundary conditions (cr£^j = erf). Hamilto¬ 
nian Eq. © is integrable at k = 0 and non-integrable for 
all K 7 ^ 0 (see e.g. BtI 148)1 1. The parameters are initialized 
to Ki, hi and then suddenly changed to K 2 , h 2 - We first il¬ 
lustrate the (double exponential) reduction in recurrence time 
for a small quench close to a non-integrable critical point. For 
small frustration k < 1/2, there is a transition from a ferro¬ 
magnetic to a paramagnetic phase increasing the external field 
h. Numerical simulations using Eq. (j^ are shown in Fig.j^left 
panel. It is evident the sharp drop close to the quantum critical 
point. 

Similar considerations also apply to the integrable case for 
which we set k = 0 in Eq. ( [T3] l. In this case the fidelity is 
precisely given by Eq. 1149 H . Again one observes a drop 


(albeit less sharp) close to the quantum critical point Fig. 
right panel. 

V. CONCLUSIONS 

We provided essentially exact formulas for the average 
recurrence time in quantum systems both for general non- 
integrable and integrable models. This is the -average- time a 
system takes to get back to the initial state up to an error 1 — u 
in fidelity. We have shown that, in the typical case, recurrence 
times are doubly exponential in V, i.e. Tr ^ (hj , 
where V is the dimensionless system’s volume normalized by 
the volume of the unit cell, J a system’s energy scale and ol 
a positive constant. Most of the time a valid approximation is 
given by Tr ^ (S/J)e“® where S is the Von Neumann’s en¬ 
tropy of the equilibrium state. For integrable systems instead 
the recurrence times are down to a simple exponential in V. 

A possibility to drastically reduce such astronomical recur¬ 
rence times is to perform a small quench experiment. By this 
we mean preparing the system in the ground state; slightly 
change some Hamiltonian parameters, and let the system 
evolve undisturbed thereafter. In this situation recurrences 
happen on a time-scale which is only exponential in the vol¬ 
ume. Furthermore if the quench is performed close to a quan¬ 
tum critical point recurrence times become roughly size in¬ 
dependent. This drastic reduction opens up the possibility to 
experimentally observe such fast recurrences in several nearly 
isolated quantum platforms. In particular ion traps with order 
of 10 atoms seems to be particularly suited. 

The author would like to thank Leonardo Banchi and Larry 
Goldstein for interesting discussions. This work was sup¬ 
ported under ARO MURI Grant No. W91 lNF-11-1-0268. 
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Appendix A: Densities of zeroes in the non-integrable case 

The fidelity is the modulus square of the characteristic func¬ 
tion 


x{t) = ^Pue (Al) 

n 

with weights pn = (tplUnlip)■ The investigation of sums of 
the form Eq. ( |A1| | was initiated long ago by Lagrange ll50l . 
The study of the variation of perihelion leads to a study of the 
variation of the argument of such a sum where the number of 
terms is the number of planets. Here instead we are interested 
in the behavior of the modulus. 

We now write x = -A -P zF such that = X"^ -P and 
P = 2XX' + 2YY' with 


X = '^pjCos{Ejt) 

j 

(A2) 

Y = '^pjSiniEjt) 

j 

(A3) 

X' = - sHEjt) 

j 

(A4) 

Y' = ^pjEjCOs{Ejt). 

(A5) 


0 


Note that, to avoid an overburdened notation, we use the 
same latter to indicate both a function of time (e.g. X{t)) 
and a random variable (X) obtained convoluting t seen as a 
random variable with uniform distribution in [0,T] and tak¬ 
ing the limit T —oo. We also use the compact notation 
X = {X,Y, X',Y') and a: = (x,y,x',y'). With the help 
of the joint distribution function Px{x), the density of zeros 
can be written as 

D{u) = J dxdydx dy S{u — x^ — y^)2 |a;a:' -P yy\ x 

Pxix,y,x' ,y'). (A6) 

From equations and the RI assumption it follows 

that X is a sum of independent random vectors. Indeed the 
characteristic function of X (fe = (^, p, rj')) is 

px{k)=e'‘^'^ (A7) 

/ d'O 

exp[zylj cos(^j) + iBj sin(79j)] (A8) 

J 

(A9) 


sin2(z9f 
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where Jq is the Bessel function of the first kind and with and the coefficients are fixed by 


Aj = pj + EjT]') and Bj = pj{rj — Ej^'). Using Eq. (A9 1 it 
may be possible to rigorously prove normality of (a properly 
rescaled version of) X under some conditions using the cen¬ 
tral limit theorem for triangular arrays. This would require to 
have a well defined family of models parametrized, for exam¬ 
ple, by the size. This is not always possible and so it will not 
be further pursued here. Discarding the small, O(fc^), term 
in Eq. (A9 1 we realize that Px{x) has approximately the fol¬ 


lowing Gaussian form 


^x/2 



tt^A ’ 

with inverse correlation matrix given by 


/ F 0 0 -E\ 

E-i - - ^ 

0 F E 0 

{DF - E^) 

0 E D 0 


\-E 0 0 D ) 


(AlO) 


having defined 


d = -£p] = ^ 

3 

3 


(All) 

(A12) 

(A13) 

(A14) 


The positivity of the matrix E reduces to A := DF — E^ > 0 
which can be easily proved using Cauchy-Schwarz inequality 
with the vectors v = (pi,p 2 , ■ ■ •) and w = {Eipi,uj 2 P 2 , ■ ■ ■) 
(under assumption of convergence) \v ■ w 1^ < ll^’f Ikf ■ 
We now pass to polar coordinates (cc, y) = 
pi (costal,sinr)i) and {x',y') = p 2 (cosi? 2 ,sin^) 2 )■ 

The density of zeroes becomes 

D{u) = j dxdydxdy' S{u — x^ — y^)2 \xx + yy\ 

X Px{x,y,x',y') (A15) 

pOO pOO 

= dpi dp2{piP2f5(u- p\)Px{p,’d) 

Jo Jo 


(A16) 


p2TZ p27T 

X / ddi / (it72 2 |cos(t?i — t) 2 ))| ■ (A17) 

Jo Jo 


Using the Gaussian approximation Eq. (AlOi, the integra¬ 
tion over Pi is trivial and the result is 


pOO n2'K ^277 

D{u) = y/u / dp2pl / ddi / dl?2 

Jo Jo Jo 

X |cOs(l?i - t?2)| .P(P2,t?l,t?2) (A 18 ) 

where 

1 


P{P2,di,d2) = e-“^2-/3sin(,?i-t)2)p2-7. 


a;E ^a:/2| ^ = 

Dpi + 2Ey/up2 sin(i7i - i? 2 ) + Fu 


DF- E"^ 


and are given by 


“= A 


3 = 


7 = 


D 

A 

2Ey/u 

A 

Fu 


A ■ 


(A19) 

(A20) 

(A21) 


Now the angular integration can be performed 


p27T p2T: 

/ ddi di ?2 |cos(t)i 
Jo Jo 

= (A22) 

3p2 

The remaining integration over p 2 is Gaussian and the final 
result is extremely simple 


D{u) = 


2 




(A23) 


Using the ensemble the variance of the en¬ 

ergy {UJ is AE^ = [EjP]E]/P - {EjP]E,/Pr] = 


/^2 


A/J^ and we obtain 




(A24) 


which, upon remembering that Tr = \/D{u), coincides with 
Eq. (§. The exponential term in Eq. (j^ could have been ob¬ 
tained with the following (wrong) argument sometimes used 
in physics. According to this argument, Tji{u), i.e., the aver¬ 
age of the solutions of P{t) = u, is approximately estimated 
by the inverse probability (density) of the event P{t) = u 
(time some unknown time-scale J~^). The probability den¬ 
sity of the fidelity has been computed in ll^l2^ and is given 
by Px(u) = t?(u) eyLp{—u/P)/P {'d{u) Heaviside function). 
This argument then predicts, Tfi{u) ^ J~^Feyip(u/F) 
which fails if u is not sufficiently away from zero. 


Appendix B: Proof of zero density for u = 0 

Throughout this section we assume that only a finite num¬ 
ber d' of pn’s are non-zero. 


-2 A 


The result ( |A23| l predicts a 
zero density of orthogonalization times, i.e. times for which 
F{t) = 0. It is legitimate to ask weather this holds exactly 
or is an artifact of the Gaussian approximation used to derive 
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Eq. ( A23| l. We then go back to equation (A15i in its exact 
form that we re-write here for clarity 

pOO n‘2'K f'‘2'K 

D{u) = y/u / dp 2 pl / ddi / |c 0 s(t?i - 1 ? 2 )| 
Jo Jo Jo 

xPxiVucosdi, ^/usm^lP 2 cost? 2 , P 2 sin 1 ^ 2 ) 

=: ^/uR{u) (Bl) 


We will show that R{u) is bounded for all u, in particular for 
u = 0. The result then simply follows from Eq. •ET) taking 
the limit m —0. Let us analyze the logic. Eirst we show 
that, by construction, Px{x) is compactly supported. Then 
we prove that Px (x) is also bounded thus implying that R{u) 
is bounded. 


Lemma: The function Px{x,y,x',y') is compactly sup¬ 
ported provided (|Tf|) is finite, more precisely, 
Px{x,y,x',y') = 0 outside the (hyper) rectangle 
|x|,|y| < land|x'|,|j/'| < C = J2jPj\Ej \ = (l-H’D- 
The latter quantity is certainly finite for the finite di¬ 
mensional case, physically it can be postulated to grow 
as the system’s volume. 


It is obvious that Px{x, y, x', y') = 0 for \x\ , \y\ > 1. On the 


other hand \X’\ 


-Y.jPjEjCOs{Ejt) 


< J2jPj \Ej\ = 


C, QED. Now Px{x) is bounded provided the characteristic 
function (its Eourier transform) Lpx{k) G The re¬ 

minder of this section is devoted to asses the conditions under 
which ipx{k) is summable. 

Remind the definition of px{k) 


d - 

^x{k) = (B2) 




where d' is the number of non-zero pj and 


/ 1 


M,=p 


0 -Ej 
\E, 0 


1 -E, 




EA 


E]} 


(B3) 


The matrix Mj is hermitian, positive semi-definite, has eigen¬ 
values 0([x2] and p|(l -f £’|)[x2] and a null space Vj = 
span 0,0, —1), (0, iJj, 1,0)}. The Bessel function is 
smooth and everywhere bounded and hence so is ipx{k). To 
see wether it is summable we must look at its behavior as 
k ^ 00 . Remember that, for x —?► 00 , neglecting an oscil¬ 
lating factor, we have | Jo(a:)| ~ ! A- Let us fix an 

e > 0 and define the cylinders (in this picture the planes be¬ 
come lines) Tj = {k\kMjk < e^}. We call “Tubes” the 
collection (union) of all the T/s. We want to see if and when 
the following integral is convergent 


Let us call R (1 2 ) the first (respectively, the second) integral 
above. Eor clarity we also go to spherical coordinates k = 

{p, n) and define ^AMjk = pfj(Q). Note that fj(yi) = 0 
on Vj but /j(n) > 0 on K'’‘\Tubes. Let us first consider /i. 


d' 


h= Y[\Jo{pfj{mdk. (B5) 

</E'^\Tubes 


To estimate convergence of the above we look at the integrand 
when p ^ 00 . We use then the asymptotic expansion of the 
Bessel function and obtain (apart from unimportant constants 
and, bounded, oscillating factors) 




1 

pd'! 2 


d' 


n 


1 

AIM' 


(B6) 


Note that the expansion in Eq. ( |B6| l is well defined because 
jj (n) is never zero in the domain of integration. Since we are 
in four dimensions we deduce that /i is convergent if d'/2 > 
4, i.e., d! > 8. Let us then look at R. By sub-additivity of the 
Lebesgue integral 


R 


iTiUTs- 

d' 

^ A 


\(fx{k)\dk < 



\pxik)\dk (B7) 

(B8) 


i=i 

Let us consider R^i, the arguments are the same for all j. On 
Ti, since AMik < , R (^\/AMih^ ~ 1 for small e, 
hence we use the simple bound | Jq (x) | < 1 and obtain 


< 


1 j=2 \ 


AMjk 


(B9) 


We also consider the change of variable k = Uf q where Ui 
is the matrix which diagonalizes Mi, i.e. UiMiUi = Di. 
Explicitly 


Ui = 




/ El 0 0 -1 \ 

0 El 1 0 

1 0 0 
V 0 1 -El 0 ) 


(BIO) 


In these variables R = {q\ + ql < /[p\{l -f Ef)]}. We 

define Mj = UiMjUi , so that 


-^2,1 < 


/Ti 


d' 


i=2 


q'^Mjq ) dq. 


(BID 


\(px{k)\dk = 

\(px{k)\dk 

R^^VTubes 


\(px{k)\dk. (B4) 


> Tubes 


Eor e non-zero, there are points q in R such that q^^Mjq — 0, 
however these regions are bounded. To see this note that the 
planes VjS intersects only at the origin. The thicked version 
of the VjS, Rs, intersect only in a bounded region around 
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the origin, because the Mj are positive semidefinite. In other 
words there is no way to make qi , (72 run to infinity while hav¬ 
ing q^MjQ = 0 such as to render the corresponding Bessel 
function ineffective. If this argument is not convincing let us 
look directly at Mjq. With a further change of variables 
q = (pcos('!?), psin(i9), i?cos((/)), i?sin(0)), with Jacobian 
Rp, we obtain 


+ 2 {El- Ej + EfEj - EiE]) Rp cos(i9 - . (B12) 

For p sufficiently large the above is never zero. In fact, for 
large p. 


^jq = ^^2 P^ + 0{p), (B13) 

such that, neglecting an oscillating factor, at leading order 


Jf 




p—»CX ) / 2 yj\ + E‘^ 1 

^ V '^P3 1^1 -^il 


whereas next to leading terms decay ever faster. All in all, ne¬ 
glecting un-important constants, we are led to study the con¬ 
vergence of the following integral (say for a certain c >0) 


|•2^T |•2■K rf-/\/pi(J+El) ncx> 2 ^ 

I "Vo 

The above is convergent for (d' — l)/2 > 2, i.e. d' > 5. 
The same arguments can be repeated for all the I 2 J imply¬ 
ing that I 2 is convergent for d' > 5. Overall is (fx{k) is in 
L^(W^) provided d' > 8. Under this condition we then have 
\Px{x)\ < f dk\(px{l^)\ < 00 . The result then follows not¬ 
ing that |i?(ti)| < C^sup^ \Px{x)\. 


Appendix C: Solutions of J^(f) = 1 

We assume here that the discrete spectrum is finite and 
J2‘j=iPj = 1- ■^(^) = 1 if and only if x(f) = e~^^, in 
other words, with E'„ = En — (j) 


d 

^p^e-**<=l. (Cl) 

Given the constraint J2j=i Pj = f > above equation can be 
satisfied if and only if all the phases are integer multiple of 27r, 
i.e. there is a time t for which E'^^t = 2'Kmn with integers 
(depending on t). Assume f ^ 0, in this case E'^jran = 
27r/f =constant. But this violates the assumption of rational 
independence. In fact, if d is even, simply take (i/2 times 
rational coefficients c„ = l/r7i„ and the other half, c„ = 
— 1 /to„. Then Yln=iPn<^n = {d/2 — d/2)2Tr/t = 0. For 


d = 2p+l odd take p times Cn = \/ m„ and p — 1 times c„ = 
— l/m„. For the remaining two cases pick c„ = —l/(2mn). 
Then Yi=i = (p - (p - 1) - 1/2 - l/2)2^/f = 0. 
But the set {E'^} is rationally independent if the set {En} is 
which is a contradiction. Hence f = 0 is the only solution of 
E{t) = 1. 


Appendix D: Density of zeroes for quasi-free fermions 

Here we consider a system of quasi-free fermions whose fi¬ 
delity is given by Eq. ([^. The first step is to find the joint dis¬ 
tribution Pz,z'{z,z') = S{Z{t) — z)S{Z'{t) — z'). Clearly 

Z'(f) = ^In [1 - afcSin^(fefc/2)] (Dl) 

k 

^ akCk sin(fefc/2) cos(fefc/2) 

^ 1 - afeSin^(fefe/ 2 ) 

which we write as Z = Zk, (Z' = z/) for ease of 

notation. We now assume rational independence of the one- 
particle energies e^. With this assumption all variables are 
independent and uncorrelated except that (z/) correlates 
with Zfe (z/). In fact it turns out that (z/) = 0 and {zkz}) = 
0 because Zk is even while z/ is an odd function of t. In 
general one expects a central limit theorem and a Gaussian 
Pz,z' {xi,X 2 )- The only non-vanishing correlations are 


In J" = ^ Zfc 

k 

k 

4' = (4)"- 


(D3) 

(D4) 

(D5) 


Time averages are translated into phase-space averages using 
the assumption of rational independence of the one particle 
energies. Evaluating the integrals one finds 


Zk = 2 In 


1 + y/l — Q-k 


(p \2 _ (2 2i/l ak Uk) 2 

V^fe/ r, n -— ^k- 

2 v 1 ~ Ctfc 


(D6) 

(D7) 


The expression for (z^)^ is slightly cumbersome. The sec¬ 
ond moment can be written as 


{Zk? 


[\n{l-aksii/^ {d/2))f ^ 
2 [ln(l-a^ 

E Jo 


(D8) 

(D9) 


where we changed variable in the second line . The last inte¬ 
gral can be expressed in terms of logarithms and dilogarithms. 





























10 


The final expression for the variance := {zkY ~ 


All in all, the joint distribution of Z, Z' is given by 


Pz,Z'{z,z') = 


2'KazCkz‘ 


exp 


{z-lnPy {z') 


l\2 


2a| 


24'J 


4 = -41og^ (Vl - afc + l) -41og(2)log(afe) 

+ 4 log (4 - 4Vl - Qffc) log (Vl - Oik + 1) 

2 \ .. . (2, [y/l — ak + l) 


4i7r log 


-4Li2 


■\/l ~ ctfc + 1 

—Oik + 2y/l — ak + 2 


4Li2 


ctfe 


(DIO) 


The evaluation of the densities of the zeroes of the equation 
In = Inu proceeds similarly as in Sec. . The evaluation 
of the integrals turns out to be simpler because the variables 
Zk and 2 ^ are factorized . The final result is 

D(u) = exp T—(Inw — In J^)^/(2cri)l , (Dll) 
-nuz 

which gives Eq. Q. 




















